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• We investigate the structures of the near-plate velocity and temperature profiles at different horizontal posi- 

tions along the conducting bottom (and top) plate of a Rayleigh-Benard convection cell, using two-dimensional 
(2D) numerical data obtained at the Rayleigh number Ra= 10 8 and the Prandtl number Pr= 4.4 of an 
Oberbeck-Boussinesq flow with constant material parameters. The results show that most of the time, and 
'■^3 ■ for both velocity and temperature, the instantaneous profiles scaled by the dynamical frame method [Q. Zhou 

and K.-Q. Xia, Phys. Rev. Lett. 104, 104301 (2010)] agree well with the classical Prandtl-Blasius laminar 
boundary layer (BL) profiles. Therefore, when averaging in the dynamical reference frames, which fluctuate 
with the respective instantaneous kinematic and thermal BL thicknesses, the obtained mean velocity and 
temperature profiles are also of Prandtl-Blasius type for nearly all horizontal positions. We further show that 
in certain situations the traditional definitions based on the time-averaged profiles can lead to unphysical BL 
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thicknesses, while the dynamical method also in such cases can provide a well-defined BL thickness for both 
£h~* the kinematic and the thermal BLs. 

A' 

I. INTRODUCTION of the turbulent RB system 3 .. Recently, the GL the- 
ory was successfully extended to the very large Rayleigh 
In a series of recent studies^ we have experimentally number regime (the so called ultimate range), in order 
and numerically analyzed the structures of the kinematic to predict the experimentally observed multiple scaling 
and thermal boundary layers (BLs) in the vicinity of the 01 the heat transfer-! 4 , and to the rotating case to pre- 
horizontal top and bottom plates, where the fluid layer dict the heat transfer enhancement. As the GL the- 
is heated from below and cooled from above in turbulent or y is based on the assumption that the BL thickness 
Rayleigh-Benard (RB) convection, in the central region scales inversely proportional to the square root of the 
of the RB cell 3 -^. The dynamics and the global features Reynolds number according to Prandtl's 1904 theory, the 
of the thermal convection system are strongly influenced, validity of Prandtl-Blasius BL flow needs to be tested 
sometimes even dominated by the properties of the BL also locally. Note that comparison of the mean bulk 
flow. Nearly all theories of the heat transport in tur- temperature calculated using the Prandtl-Blasius thc- 
V bulent RB convection, from the early marginal stability 01 T with that measured in both liquid and gaseous non- 
theory^ to the Shraiman & Siggia (SS) model&Z and to Oberbeck-Boussinesq RB convection shows very good 
P. ! the Grossmann & Lohse (GL) theory 8 .--!!, are essentially agreement i 8 .. In addition, the kinematic BL thick- 
BL theories. Therefore, it is a key issue of turbulent RB ness evaluated by solving the laminar Prandtl-Blasius BL 
convection, how the near-plate velocity and temperature equations was found to agree well with that obtained in 
profiles look like the direct numerical simulation (DNS)!^. 

Specifically, the GL theory has achieved great success 
in predicting the global quantities, such as the Rayleigh 

number— and Prandtl number-H dependence of the heat A. BLs along the cell's central axis 

flux, i.e., the Nusselt number, and amplitude of the large 

scale circulation (LSC), i. e., the Reynolds number, Previous works about the kinematic and thermal 

BLs!* 2 - mainly focused on the cell's central vertical axis. 
Although the flow in the bulk is turbulent, the BLs are 
considered to behave still laminar at least scaling wise 
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(see Fig. [T^b) of this paper for the values of the shear 
Reynolds numbers in the present case). Indeed, in a 
time-averaged sense, it was found experimentally that the 
kinematic BL thicknesses A„ near the sidewalls of a cubic 
cell^S and near the bottom plate of a rectangular cell 21 
obey the Prandtl scaling for a laminar flat plate BL, i.e. 
X v ~ -Re -1 / 2 , where Re is the Reynolds number of the 
LSC in the RB system. Furthermore, certain wall quan- 
tities, such as the wall shear stress, the friction velocity, 
and the viscous sublayer thickness, were also found to fol- 
low the Prandtl scaling 21 .. However, direct comparisons 
of experimental velocity^ and numerical temperature 2 ^ 
profiles with the respective classical Prandtl-Blasius pro- 
files show significant deviations, especially for the dis- 
tances from the plate around the BL thickness. It was 
argued that such deviations should be attributed to the 
intermittent emissions of thermal plumes from the BLs 
and the corresponding temporal dynamics of the BL a 1 ' 24 . 
This led to the study of the BL structures in dynam- 
ical reference frames, which fluctuate with the instan- 
taneous BL thicknesses, rather than in the laboratory 
framed. When resampling the velocity and temperature 
fields' data in such dynamical frames, both the mean ve- 
locity and the temperature profiles were found to agree 
well with the respective theoretical Prandtl-Blasius lam- 
inar BL profiles over a wider parameter range of both 
Ra and Pr-i^. Moreover, when the instantaneous veloc- 
ity and temperature profiles are rcscalcd by their respec- 
tive instantaneous BL thicknesses, it was found that the 
Prandtl-Blasius profiles not only hold in a time-averaged 
sense, but are most of the time also valid in an instanta- 
neous sensed. A dynamical BL rescaling method has thus 
been established, which extends the time-independent 
Prandtl-Blasius BL theory to the time-dependent case, 
in the sense that it holds locally at every instant in the 
frame that fluctuates with the local instantaneous BL 
thickness. All this was, as mentioned, shown for the cen- 
ter range of the RB cell. 
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FIG. 1. (color online). The time-averaged vector map of the 
whole velocity field v = (u,w) (see Sec. Ill B I for the details). 
For clarity, a coarse-grained vector map of size 26 x 50 mesh- 
points is shown. The magnitude of the velocity v = \/u 2 + w 2 
is coded in both color and the length of the arrows in units 
of m/s. The time average is taken over a period of 80 min 
corresponding to 30 000 velocity maps and to 500 large eddy 
turnovers (LET), with Tlet = (2H + 2D)/u LS c = 9.6 s, 
ulsc = 0.0I7ms _1 , and 60 velocity maps per LET. If one fol- 
lows the stream trace passing through the maximal velocity 
ulsc and uses the numerically measured local velocities along 
the approximately elliptically shaped circumference of the 
stream strace, the turnover time is only 7.2 s, corresponding 
to 667 LET in 80 min. - The solid curve marks the kinematic 
BL thickness near the bottom plate. The Reynolds number 
of the LSC is ReLSC = ulscH/v = 1036, where ulsc is the 
maximal velocity magnitude of the LSC, and the Reynolds 
number of the lower left corner roll is Re cr = u cr i cr /v — 134, 
where u cr = 0.01 1 ms~ is the maximal velocity magnitude of 
the lower left corner roll and £ cr ~ 0.2H is the typical length 
scale for the corner roll. 



B. The spatial structures of the BLs 

As a closed system, turbulent thermal convection in 
an RB cell develops rather complicated flow structures, 
partly due to the interactions between the flow and the 
solid walls, cf^. The kinematic and thermal BLs along 
the cell's central vertical axis thus cannot reveal all BL 
properties, especially not for the BLs in the regions near 
the cell's corners. In Fig. [TJ we show an example of 
the time-averaged vector map of the whole velocity field 
obtained from a two-dimensional (2D) simulation with 
the Rayleigh number Ra= 10 8 and the Prandtl number 
Pr= 4.4 (for the details of the simulations we refer to 
Sec. HI]) . As usual the Rayleigh number is defined as 
Ra = agH^/^/vK and the Prandtl number as Pr= u/k. 
Here v, k, a, and g are the kinematic viscosity, thermal 
diffusivity, isobaric thermal expansion coefficient, and 
gravitational acceleration; H denotes the height of the 



container and A the temperature difference between the 
hotter bottom and the cooler top temperature. It is seen 
clearly that the overall flow pattern is an counter- clock- 
wise rotatory motion. While in a three-dimensional (3D) 
cylindrical cell2£ the mean flow was found to be ellip- 
tically shaped, the large-scale circulation (LSC) in the 
present case looks a bit more stadium-like shaped with 
its long and short axes pointing approximately to the 
cell's two diagonals. There are several smaller secondary 
rolls at the four corners of the cell: two larger clock- 
wise rolls at the two opposite corners adjacent to the 
short axis of the LSC ellipse and much smaller vortices 
at the two opposite corners adjacent to the long axis of 
the LSC ellipse. Thus the flow near the horizontal plates 
can be divided into the two corner-roll regions and the 
central region dominated by the LSC. To see this more 
clearly, we plot in Fig. [5] the horizontal profile of the 
time- averaged horizontal velocity u(x) near the bottom 
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FIG. 2. (color online). The time-averaged horizontal velocity 
profile u(x) as a function of x/D obtained near the bottom 
plate (z/H = 0.0036). The profile can be divided into three 
regions: the left corner roll (region I), the LSC (region II), 
and the right corner rolls (region III). The vertical solid lines 
mark the boundaries between the three regions at x a and xj,- 
Note that at Ra = 1 x 10 s both the kinematic and thermal BL 
thicknesses are larger than 0.009// for all horizontal positions, 
see the solid curve in Fig. [7] (a) . Thus the horizontal velocity 
profile in this figure is well within both BLs. 



i.e., Xth tends to eventually become uniform along the 
plate at very large Ra. This behavior can be shown to 
result from the shape evolution of the LSC. Namely, its 
shape evolves from a tilted and nearly elliptical shape at 
low Ra to a more squarish shape at high R a 25 ' 28 . The 
squarish-shaped LSC at high Ra will make the mean flow 
near the horizontal plates to be more parallel to the plates 
and hence leads to more uniformity of the BLs. 

The spatial structures of the kinematic BL in water 
was first studied experimentally by Qiu and Xia3£ in 
a cubic cell. It was found that the magnitudes of the 
LSC, the shear rate, and the kinematic BL thickness all 
change significantly across the horizontal plates both par- 
allel and perpendicular to the LSC. Direct comparison 
between the observed temperature BL profiles and the 
Prandtl-Blasius thermal profiles at different horizontal 
positions were performed by Sugiyama et al^ in 2D and 
by Stevens et al£Q in 3D numerical simulations. It was 
found that due to the rising (falling) plumes near the 
sidcwalls the deviations of the numerically calculated BL 
profiles from the Prandtl-Blasius profiles increase from 
the center of the horizontal plates towards the sidewalls. 

C. The objective of the present work 



plate (z = 0.0036-ff). One can distinguish three different 
ranges of x that differ from each other by different val- 
ues of u(x). Region I is the left-corner- roll region where 
the flow is dominated by the clockwise corner roll and 
u(x) is negative, region II is the central region where 
the flow is dominated by the counter-clockwise LSC and 
u(x) is positive. The transition point between region I 
and II is identified as x a /D = 0.44. In region III the flow 
is somewhat complicated and both negative and positive 
u{x) are observed, indicating several small vortices in the 
right corner. The transition point between region II and 
III is identified as x^/D = 0.77. - On the upper plate 
the velocity profile is correspondingly, but from right to 
left. 

Such complicated flow structures near the horizontal 
plates highlight the need to study the horizontal depen- 
dence of the local BL profiles, both of the velocity and the 
temperature. Compared to the large amount of studies 
on the near-plate velocity and temperature profiles along 
the cell's central vertical axis, however, studies on the 
spatial dependence of these profiles off the center line are 
very limited both experimentally or in simulations. 

The spatial structure of the thermal BL in water in 
the range 10 8 < Ra < 10 10 has systematically been stud- 
ied first by Lui and Xia in a cylindrical cell2£ and then 
by Wang and Xia in a cubic cell^l. Both experiments 
have shown that the thermal BL thickness above the 
bottom plate, Xth, depends on the horizontal position x 
along the plate, and the scaling exponent of Xth with Ra 
varies between —0.35 and —0.28. However, this position- 
dependence is expected to decrease with increasing Ra, 



In this paper we want to extend previous worksi^ deal- 
ing with the plate's center to the whole bottom (top) 
plate with the help of 2D DNS. Our results will show that 
the idea of the dynamical BL thickness rescaling method 
works well for almost all horizontal positions, i.e., the 
mean BL profiles obtained at nearly all horizontal posi- 
tions can be brought into coincidence with the Prandtl- 
Blasius laminar BL profiles, if they are re-sampled in the 
time-dependent frames of the local BL thicknesses, for 
both velocity and temperature. 

II. DEFINITIONS, NUMERICAL PARAMETERS, AND 
DATA ANALYSIS 

A. Numerical methods 

The mathematical model, the numerical scheme, and 
the code validation have been described elsewhere^. 
Thus we give only their main features here. The compu- 
tational domain consists of a 2D square cell of horizontal 
length D = 4.078 cm and hight H = 4.078 cm, the aspect 
ratio is thus T = D/H = 1. The flow is calculated by 
numerical integration of the 2D time-dependent imcom- 
pressible Obcrbeck-Boussincsq equations with a fourth- 
order finite-difference scheme. No-slip velocity boundary 
conditions are applied to all four solid walls. As temper- 
ature boundary conditions the two sidewalls are chosen 
to be adiabatic (no flux) , while at the colder top and the 
warmer bottom plates the temperatures are fixed. The 
mean temperature is chosen as T m = 40° C and water as 
the working fluid. Then the kinematic viscosity, thermal 
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FIG. 3. (color online). Time trace of Q(t). 



diffusivity, and isobaric thermal expansion coefficient are 
v = 6.6945 x 1CT 7 nrV 1 , re = 1.5223 x 1CT 7 naV" 1 , 
and a = 3.8343 x 10~ 4 K _1 . The resulting Prandtl num- 
ber is Pr = v/k = 4.4. During the computation, the 
temperature difference across the fluid layer was fixed at 
A = 40 K. The corresponding Raylcigh number then is 
Ra = agH 3 A/vn = 10 8 . 

For these values of the control parameters the Nusselt 
number Nu (in its usual definition as a z-independent 
area average; here in DNS after z-averaging in addition) 
is Nu = 25.62, which corresponds to the x-independent 
overall thermal BL thickness X th /H = l/2Nu = 0.0195. 
If Nu is calculated with the slopes of the area (here the 
x) averages d z (6) x at the bottom or top plates, the re- 
spective values are 25.65 and 25.69, being within 0.3% 
with the z-average. 

We denote the computational domain as the (x, z)- 
plane. Then at any time t the horizontal and vertical ve- 
locity components, u(x, z, t) and w(x, z,t), and the tem- 
perature 9(x, z, t) are obtained. For the temperature field 
a non-dimensional temperature Q(x, z,t) is introduced, 
namely, 



Q(x, z, t) 



e 



hot 



(x,z,t) 



A/2 



(1) 



Here Obot is the fixed hot temperature at the bottom 
plate. Based on this definition, the fixed dimension- 
less temperatures of the top and bottom plates are 
Q(x,H,t) = 2 and Q(x, 0, t) = 0, respectively, and the 
mean bulk temperature is O = 1. 



B. Reversals of the LSC 

Previous studies have shown that for the numerical 
parameters used in the present work, the mean flow of 
the 2D RB convection would experience spontaneous flow 
reversals, due to the competitions between the corner 
rolls and the LSC—. This can be characterized by the 
global angular velocity of the whole flow field, defined as 



w \x~0.5D 



u(x, z, t) \ 
Z-0.5H/, 



(2) 



where (■■•)« denotes the spatial, i.e., (x, z) average. 
Based on the definition of fl(t), Q(t) > indicates the 
counter-clockwise rotation of the mean flow while the 
mean flow rotates clockwise when ft(t) < 0. The events 
of flow reversals can therefore be identified through a 
sign change of Cl(t). Figure [3] shows a time series of f2(t). 
One can identify a total of 32 reversals in a period of 
80 minutes. Such a flow with so many reversals in the 
relevant averaging time yields a nearly zero time mean ve- 
locity field, because the counter-clockwise and clockwise 
velocity contributions cancel under time-averaging. To 
overcome the influence of the reversals, before analyzing 
the BL data we transform the 2D velocity and temper- 
ature fields, from u(x, z, i), w(x,z,t), and Q(x,z,t) to 
u'(x,z,t), w'(x,z,t), and 0'(x, z,t), as follows: 

1. at any time t, if O(t) > (counter-clockwise), the 
2D velocity and temperature fields are taken un- 
changed, i.e. u'(x,z,t) = u(x,z,t), w'{x,z,t) = 
w(x, z, t), and Q'(x, z, t) = 0(x, z, i); 

2. at any time t, if f2(t) < (clockwise), the 2D 
velocity and temperature fields are reflected with 
respect to the axis x/D = 0.5, i.e., u'(x, z,t) = 
—u(D — x,z,t), w'{x,z,t) = w(D — x,z,t), and 
Q'(x,z,t) = e(D-x,z,t). 

After such transformation, the LSC (if it exists; dur- 
ing the reversals the LSC breaks down—) would ro- 
tate counter-clockwise for all times t. A resulting time- 
averaged velocity field is illustrated in Fig. Q] For 
convenience, in the remainder of the paper we also use 
the notations of u(x, z,t), w(x,z,t), and Q(x, z, t) for 
the velocity and temperature after the described trans- 
formation. Since the governing equations are strictly 
Oberbeck-Boussinesq (temperature independent fluid pa- 
rameters), one expects top-bottom symmetry. Though 
for theoretical reasons this must be valid, we also have 
checked that numerically for the velocity field explicitly. 
We found complete agreement in all details of the in- 
dependent structures. All respective profiles collapse un- 
der top-bottom and x D — x mapping. We thus need 
to discuss only the velocity and the temperature profiles 
near the bottom plate (without labeling that, since it also 
holds for the top plate, respectively). 



C. The kinematic and thermal BL thicknesses obtained 
from the time-averaged profiles 

With the measured time-averaged velocity and temper- 
ature fields, the kinematic and thermal BL thicknesses 
X v (x) and Xth{x) can be defined, respectively, via the 
z-slopes of the time-averaged velocity and temperature 
z-profilcs for each horizontal position x, 

u(x, z) — (u(x, z, t)) and Q{x, z) = (0(x, z, t)}, (3) 

cb£i22, where (■ ■ • ) denotes the time average. Usually the 
width A^J 9% is considered instead. Since in RB geometry 
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FIG. 4. (color online), (a) Normalized kinematic and thermal 
BL thicknesses in z-direction near the bottom plate, X v (x)/H, 
\™ ax (x)/H, and \ t h(x)/H, as functions of the horizontal po- 
sition x, obtained from the time-averaged velocity and tem- 
perature fields u(x,z) and Q(x,z). The vertical solid lines 
mark the boundaries x a and Xb between the three regions. 
Note the different scales used for the kinematic and ther- 
mal BL thicknesses, (b) The time-averaged velocity z-profile 
u(x,z) for fixed x/D = 0.447 (near the boundary x a ). The 
tilted dashed line is a linear fit to the near-plate part of the 
it-velocity's z-profile and the horizontal dashed line marks the 
maximum horizontal velocity. 



the velocity does not level asymptotically but decreases 
again with z, we here consider \™ ax or XYjf x , defined 
as the width at which the profile reaches its (first) maxi- 
mum. The max- and the 99%- widths differ only minutely. 

Figure QJa) shows X v (x)/H and Xth{x)/H as functions 
of x/D. The thermal BL is much thicker at the two 
corners because of the rising plumes near the sidewalls. 
The horizontal distribution of Xth{x) is asymmetric and 
there exists a minimum value of Xth{x) at x/D ~ 0.4. 
These features of the thermal BL thickness as a function 
of x are similar to those observed in a cubic cell^I and to 
those in a numerical study of a square cell^. In contrast, 
in a cylindrical cell a symmetric "V" shape of the Xth (x) 
profile was found^. This difference was attributed to the 
effects of the sharp corners in the cubic or square cells on 
the flow2£, suggesting an impact of the cell geometry on 
the characteristics of thermal BL structures in turbulent 
RB convection. 

The x-dependence of A„ (x) is much more complicated. 



Here one finds large discontinuities and even gaps in the 
measured X v (x) profile around the two boundaries x a and 
Xb of the three regions. These features are caused by the 
strong competition between the corner rolls and the LSC 
at these positions and apparently are unphysical. The 
fluctuating instantaneous border between the corner rolls 
and the LSC makes the velocity change its sign [see Fig. 
(5Jc)] and slope near the plate very quickly, so that a kine- 
matic BL thickness cannot be well defined. Figure BJb) 
shows an example of the time-averaged velocity z-profile 
obtained at fixed x/D = 0.447 (near the boundary x a )- 
At this x-position, the strong competition between the 
corner and center rolls makes the near-plate velocities 
sometimes positive (when the flow is dominated by the 
LSC) and sometimes negative (when the flow is domi- 
nated by the left corner roll). The positive and negative 
velocities cancel each other when calculating the time- 
averaged velocities, yielding a rather small velocity slope 
near the plate (see the tilted dashed line in the figure) and 
thus a very large value of X v (x), which is obviously un- 
physical. Therefore, the traditional definition of the kine- 
matic BL thickness from the slope of the time-averaged 
z-profile cannot properly handle the situation as shown 
in Fig. |4j One may define, instead, the kinematic BL 
thickness as the distance X™ ax (x) to the near-plate ex- 
tremal horizontal velocity. However, still the sign change 
of the velocity along the horizontal direction [see Fig. [5] 
produces unphysical discontinuities and large jumps in 
the X™ ax (x)/H profile [see Fig. Ufa)]. As we shall see 
below, these unphysical results for the time-averaged BL 
thickness can be eliminated if we use the instantaneous 
BL thickness based on the instantaneous velocity profiles. 

We have also checked the profiles and thicknesses of the 
rms-temperature and -velocity profiles (cf<21). X™ s {x) 
and the slope thickness Xth (x) show the same features as 
functions of x, the rms-thickness being about 20% below 
the slope thickness. The velocity field's BL thickness 
Xl ms (x) is much closer to X™ ax (x), within 5%; its in- 
dependence reflects the same structures, including the 
discontinuities. Only the boundaries x a ,b are somewhat 
shifted. We thus do not consider the rms-fields further. 



D. The instantaneous kinematic and thermal BL 
thicknesses 

Figure EJa) shows examples of the instantaneous hor- 
izontal velocity's u(x,z,t) vertical profiles versus z for 
fixed horizontal position x and time t, in particular at 
x/D = 0.2 and x/D = 0.6. As shown in Figs. Q] 
and [2] in the time-averaged sense position x/D = 0.6 
belongs to region II where the flow is dominated by 
the LSC. The properties of u{x = 0.6D,z,t) are simi- 
lar to those reported in our previous works^. Namely, 
u(x = 0.61?, z, t) increases very quickly from to the in- 
stantaneous maximum velocity within a very thin layer 
above the bottom plate and then decreases slowly in the 
bulk region of the closed convection cell. In contrast, the 
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FIG. 5. (color online), (a) Examples of the instantaneous 
horizontal velocity's u(x,z,t) z-profiles obtained at x = 0.2D 
(red circles) and x = 0.6D (blue triangles), (b) An enlarged 
part of the velocities' z-profiles near the bottom plate. The 
tilted dashed lines are linear fits to the linear part of the veloc- 
ity profiles near the bottom plate and the horizontal dashed 
lines mark the instantaneous minimum (for x = 0.2D) or max- 
imum (for x = 0.6D) horizontal velocities near the bottom 
plate. The distances of the crossing points from the plate de- 
fine the instantaneous local kinematic BL thicknesses 5 v (x,t), 
either of the LSC or of the corner roll. 
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FIG. 6. (color online), (a) Examples of the normalized in- 
stantaneous local temperature z-profiles Q(x,z,t) obtained 
at x = 0.2D (red circles) and x — 0.6D (blue triangles), (b) 
An enlarged part of the temperature profiles near the bottom 
plate. The tilted dashed lines are linear fits to the linear parts 
of the temperature z-profiles near the bottom plate and the 
horizontal dashed line mark the bulk temperature = 1. The 
distances of the crossing points from the plate define the in- 
stantaneous local thermal BL thicknesses Sth(x,t). The latter 
one here is smaller for x = 0.2D than for x = 0.6D. 



,z-profile u{x = 0.2D, z,t) in region I shows quite different 
features. Here the flow is dominated by the left corner 
roll. It is seen that u(x, z, t) as a function of z first drops 
very quickly from to the instantaneous minimum ve- 
locity near the bottom plate, i.e., the velocity very near 
to the plate is negative, the flow is towards the corner. 
After reaching its near-plate (negative) minimum value, 
u(x = 0.2_D,z,i) rises very quickly from negative values 
through to the instantaneous maximum velocity and 
then decreases slowly in the bulk region. The kinematic 
BL at this position is produced and stabilized by the vis- 
cous shear of the corner roll, rather than by the LSC. 
The total thickness of the corner or secondary roll at this 
time is of the order 0.477. Figure [5jb) shows the en- 
larged near-plate parts of the two velocity profiles. One 
sees that both profiles are linear near the bottom plate 
and thus the instantaneous local kinematic BL thickness 
can be defined as the distance from the plate at which 
the extrapolation of the linear part of u(x, z, t) versus z 
crosses the horizontal line passing through the instanta- 
neous near-plate cxtrcmum horizontal velocity (e.g. min- 



imum for x = 0.2D and maximum for x = 0.6D). We 
denote these instantaneous local BL thicknesses by the 
symbol S v (x,t) to distinguish them from the thicknesses 
X v (or X t h) of the time-averaged z-profiles. The dashed 
lines in Fig. [5{b) illustrate how to determine 5 v (x,t) as 
the crossing point distances for the two ^-positions. 

Figure^a) shows examples of the z-profiles of the nor- 
malized instantaneous local temperature 0(x, z,t), ob- 
tained at fixed x = 0.2D and x = 0.6D and chosen time 
t. Unlike the case of u(x, z, t) in Fig. [5l although the 
two positions are dominated by different flow directions, 
both 0(x, z,t) here exhibit similar features as functions 
of z. &(x, z, t) increases very quickly within a very thin 
layer above the bottom plate and stays nearly constant 
at the mean bulk temperature = 1, valid in the bulk 
region. Figure [BJb) shows an enlarged near-plate part 
of the two temperature profiles. Linear pieces can be 
seen in both profiles near the bottom plate. The instan- 
taneous local thermal BL thickness (again denoted by 
the symbol 5 instead of A) 5 t h{x,t) can then be defined 
as that distance from the plate, where the extrapolation 
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FIG. 7. (color online), (a) The horizontal (independence of 
the time-averaged local instantaneous BL widths (S v (x,t)) /H 
and {5 t h(x,t)) / H . The dashed curves mark the horizontal 
dependence of Xth(x)/H and X v (x)/H for comparison, (b) 
The time-averaged values of the instantaneous shear Reynolds 
number (Re s (x,t)) as a function of x/D. The vertical solid 
lines mark the boundaries x a and Xb between the three regions 
I, II, III. 



of the linear part of Q(x,z,t) crosses the horizontal line 
passing through the mean bulk temperature, even if the 
measured temperature stays below that, as for x = 0.2D. 
The dashed lines in Fig. |UJb) illustrate how to determine 
5th{x,t) as the crossing point distances. 

As discussed in Sec. Ill C[ the kinematic local BL thick- 
nesses, X v (x) and \™ ax (x), based on the time-averaged 
velocity profiles have some unphysical features. Specif- 
ically, Fig. HJb) illustrated the limits of the traditional 
definition of the BL thickness. We now show that the 
use of the instantaneous BL can avoid some problems. 
We consider the time-averaged mean values of the in- 
stantaneous local BL thicknesses (5 v (x, t)) and (5th(x, t)) 
as typical measures of the local kinematic and thermal 
BL thicknesses, respectively. Figure [7] shows (S v (x,t)) /H 
and (<5th,(a;, t)) / H as functions of x/D. For compari- 
son, X t h(x)/H and X v (x)/H obtained from the time- 
averaged temperature and velocity profiles are also plot- 



ted as dashed curves in the figure. For temperature, one 
sees that the horizontal dependences of (5 t h(x, t))/H and 
Xth(x)/H share the same trend, with (Sth{x,t)) /H only 
a little bit larger than Xth{x)/H . For velocity, the situa- 
tion is quite different. Note that the time-averaged local 
instantaneous BL width (5 v (x, t)) / H varies smoothly be- 
tween 0.008 and 0.04 along the whole horizontal plate, 
especially around the boundaries x a and Xb between the 
different regions. No gaps exist and no extremely large 
thicknesses. This is because at any instant the instanta- 
neous local kinematic BL thickness S v (x,t) can be well 
defined regardless whether the local flow is dominated 
by the LSC or by the corner rolls at that particular in- 
stant. Therefore, (5 v (x,t)) / H can be used to charac- 
terize the typical length scale of the kinematic BLs in 
situations where X v (x)/H and X™ ax (x)/H, based on the 
time-averaged velocity profiles, are no longer capable of 
producing physically meaningful results. 

Let us briefly compare with the conventional BL thick- 
ness. The area- (here the cc-)averagcd thermal BL thick- 
ness is (Sth(x)) x — 0.024£f; this is slightly but definitely 
different from the global thickness X t h = H/2Nu = 
0.0195iJ given in Sec. II. A. The x-averaged rms-profile 
thickness (Aj™ s (x)) 2; = 0.0176iJ turns out to be less. 
- For the kinematic BL, in contrast, due to the corner 
rolls and the corresponding back flows, an x-averaged v- 
profilc is not meaningful. But one might x-averagc the 
local thicknesses to obtain (5 v (x)) x = 0.02Q7H. Com- 
pare this average with the x-dependent thicknesses in 
Fig. 4(a) and find that under averaging the smaller cor- 
ner roll BLs reduce the center roll BL thickness consider- 
ably. - Note that the slope thickness of the center roll is 
of order A t , « a/V Rs-lsc with a ~ 1.6. This differs from 
the value a = 0.5 valid in the Prandtl law for the global 
Re number, cf£~— . 

The instantaneous local BL thicknesses also yield the 
so-called local shear Reynolds number, based on the local 
kinematic BL thickness as the characteristic length scale. 
We define it as 



where u(x,S v ,t) is the instantaneous local velocity at 
z = 5 v (x, t). Figurc[7]Jb) shows the horizontal variation of 
the time-averaged values of the instantaneous local shear 
Reynolds numbers, (Re s (x,t)) . The x-profile of this time 
averaged local shear Reynolds numbers shows two peaks. 
One peak is within the lower left corner roll and the other 
is within the LSC. These local shear Reynolds numbers 
(Re s (x, t)) are for all horizontal positions x much smaller 
than the critical value Re s = 420 for the instability of 
the boundary layer, which had been proposed in the lit- 
erature cL 3 -^. This suggests a still laminar though tem- 
porally fluctuating BL for the present control parameters 
Ra and Pr. Note that in a previous experimental study^i 
it was found from the extrapolation of the Re s versus Ra 
scaling that for Pr = 4.3 a turbulent BL is expected 
to occur at Ra~ 2 x 10 13 (the similar conclusion was 
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also made in an experimental study of local heat flux 
measurements^) . 



E. The dynamical BL rescaling method 

With the measured 5 v (x,t) and Sth(x,t), the local dy- 
namical kinematic and thermal BL frames at different 
horizontal positions x along the bottom or top plates 
can now be constructed. We define the time-dependent 
relative vertical distances z*(x,t) and z^ h {x,t) from the 
plate with respect to 6 v (x,t) and 5th{x,t), respectively, 
as 



z*(x,t) 



S v (x,t) 



and Zt h (x,t) 



Sth(x,t) 



(5) 



The mean local velocity and temperature profiles, 
u*(x,z*) and Q(x,z^ h ), in the respective dynamical BL 
frames at any horizontal position x are then defined as 



u*(x,z*) = (\u(x,z = z*6 v (x,t),t)\) 



and 



9*(av 



■lh) 



(Q(x,z = Zt h S t h(x,t),t)), 



(6) 



(7) 



i.e. time-averaging over all values of \u(x,z,t)\ and 
6(2;, z,t) that were measured at different discrete times 
t, but at the same rescaled positions z* and z^ h , respec- 
tively. Here we use the absolute values \u(x, z, t)\ when 
calculating the mean velocity profiles, because u(x, z, t) 
has different signs at different horizontal positions x and 
at different time t, especially for the positions around x a 
and Xb, and u(x,z,t) with different signs would cancel 
each other partially. 

In order to characterize the shapes of the (time- 
averaged or instantaneous) local velocity and tempera- 
ture z-profiles and to study their agreement or deviations 
from the respective Prandtl-Blasius profiles, both quan- 
titatively, we compute the local shape factors Hi{x) of 
the profiles, cfi^S 



H i (x) = Sf(x)/SY l {x), with i = v,th, 



(8) 



where Sf(x) and Sj n (x) are the local displacement and 
momentum thicknesses of the profiles, respectively, de- 
fined as, 



Y(x,z) 
[Y(x,z)] m , 



dz 



(9) 



and 

TO 



Y(x,z) 



[Y(x,z)} r 



Y(x,z) 
[Y(x, z)]m, 



dz. 
(10) 

Here, Y{x,z) = \u(x, z)\ is horizontal local velocity's z- 
profilc if i = v, and Y(x, z) = 0(x, z) is the correspond- 
ing local temperature's z-profilc if i — th. As suggested 



in our previous studies^, all z-integrations are evalu- 
ated only over the range from z — to the position of 
the first maximum of the z-profiles instead of oo. The 
shape factor of a profile describes how fast the profile 
approaches its asymptotic value. The larger the shape 
factor is, the faster the profile runs to its asymptotic 
level. The shape factor of the Prandtl-Blasius velocity 
profile is H BB =2.59 independent of Pr and that of the 
Prandtl-Blasius temperature profile is H^ B = 2.61 for 
the present Pr= 4.4. The deviations of the numerical z- 
profiles from the respective Prandtl-Blasius profiles can 
then be measured by 



PB 



6Hi = Hi - H 



III. VELOCITY PROFILES NEAR THE PLATE 



(11) 



To reveal the horizontal (x-)dependence of the BL 
structures using the dynamical rescaling method, we fo- 
cus on nine different horizontal positions at x/D = ixO.l 
with i = 1, 2, and 9. As shown in Figs. Q] and [21 in 
the time-averaged sense the positions x/D = 0.1, 0.2, 0.3, 
and 0.4 belong to region I, where the flow is dominated 
by the left corner roll. The positions x/D = 0.5, 0.6, 
and 0.7 belong to region II, where the flow is dominated 
by the LSC. The positions x/D = 0.8 and 0.9 belong to 
region III, where the flow is dominated by the small right 
corner rolls. To reduce the data scatter, before applying 
the dynamical rescaling method we coarse grain the local 
horizontal velocity and temperature profiles u(x, z, t) and 
0(x, z, t) obtained at each position i and at each discrete 
time t by averaging them along the a;-direction (horizon- 
tally) over the range ix 0.1 -0.01 < x/D < ix 0.1 + 0.01. 

Figures [5Ja) and (b) show the magnitudes of the z- 
profilcs of the time-averaged local horizontal velocity 
\u(x, z)\, obtained respectively in regions I and II. Here, 
the velocities are normalized by their respective near- 
plate maximum values [|u(x, z)|] moa; and the vertical dis- 
tance z is normalized by the local kinematic BL thickness 
X v (x). For comparison, we also plot the Prandtl-Blasius 
velocity z-profile (the dashed lines), the initial slope of 
which is matched to those of the measured profiles. It is 
seen that the time-averaged velocity z-profiles measured 
at positions within region II deviate significantly from 
the Prandtl-Blasius profile (Fig. EJb)). On the other 
hand the z-profiles in region I, where the flow is domi- 
nated by the left corner roll, match the Prandtl-Blasius 
one much better (Fig. Ufa)). This observation was some- 
what to our surprise, since for the temperature profiles 
in both 2D21 and simulations increasing deviations 
from the Prandtl-Blasius profile were found when away 
from the cell center, due to the rising plumes close to 
the sidewalk Our present result may be understood from 
Fig. [TJ in which one sees that the mean flow (in the cor- 
ner roll) in region I is essentially parallel to the plate, 
which is just the case treated by the Prandtl-Blasius BL 
theory, i.e., a horizontal flow over a flat plate. In con- 
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FIG. 8. (color online), (a, b) The absolute values of the 
time-averaged local horizontal velocity z-profiles, \u(x, z)\, as 
functions of the normalized distance z/X v (x) obtained at (a) 
x/D = 0.1, 0.2, 0.3, and 0.4 (region I) and (b) x/D = 0.5, 
0.6, and 0.7 (region II). Here, \u(x, z)\ is normalized by its re- 
spective maximum value near the bottom plate, [\u(x, z)\] m ax- 
The dashed lines indicate the Prandtl-Blasius velocity profile 
for comparison, (c) The time-averaged local horizontal ve- 
locity 2-profiles u(x,z), normalized by the respective maxi- 
mum horizontal velocity near the bottom plate [u(x, z)] max , 
as functions of z/H obtained at x/D = 0.8 and 0.9 (region 
III). 



trast, the flow enters region II from above not parallel to 
the plate and only becoming more horizontal afterwards. 

Figure EJc) shows the time- averaged local velocity pro- 
file u{x, z) including its sign, measured at x/D = 0.8 and 
0.9, i.e., in region III. The velocities u(x, z) are still nor- 
malized by their respective maximum values, while z is 



FIG. 9. (color online). PDFs of the shape factors H v (x,t) 
of the rescaled instantaneous velocity profiles, see figure If 01 
obtained at (a) x/D = 0.1, 0.2, 0.3, and 0.4 (region I), (b) 
x/D = 0.5, 0.6, and 0.7 (region II), and (c) x/D = 0.8 and 
0.9 (region III). The dashed lines mark the shape factor of 
the Prandtl-Blasius velocity profile for comparison. 



normalized here by the cell's height H instead of the lo- 
cal BL widths. Due to the very small right corner rolls 
in region III, the two profiles first drop a little from 
and then increase to their maximum values. We find 
that there are respectively 1 and 4 data points between 
and the near-plate minimum velocities of the two pro- 
files. These numbers of data points are so small that it 
is meaningless to define BL thicknesses for these profiles. 
In the present study, the kinematic BL thickness was cal- 
culated only, when the number of data points between 
and the near-plate extremal velocities was larger than 5. 
Having analyzed the profiles of the time-averaged lo- 
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FIG. 10. (color online). Comparison between velocity profiles obtained at x/D = 0.1, 0.2, 0.9 near the bottom plate: 
dynamical u*(z*) (red solid lines), laboratory u(z) (blue solid lines), and the Prandtl-Blasius laminar velocity profile (black 
dashed lines). Note that the Prandtl-Blasius profile per construction stays constant once it has reached its asymptotic value, 
while the RB flow profiles decrease towards the bulk. Thus agreement can only be expected in the very BL region. As the plots 
show, in the BL range the dynamically rescaled instantaneous local profiles are very well consistent with the Prandtl-Blasius 
shapes. 



cal velocities with the tools of the time-averaged fields 
and have noted considerable deviations from the Prandtl- 
Blasius behavior, we now consider its comparison with 
the z-profiles of the instantaneous local horizontal veloc- 
ity. At first we evaluate the local instantaneous shape 
factors H v (x,t) to characterize the respective z-profiles 
of the horizontal velocity near the plate, see Fig. [9] 
Then the dynamically rescaled z-profilcs themselves will 
be considered, see Fig. [TU] The corresponding analysis 
for the thermal profiles will be presented in Sec. IV. 

Figure IHJa) shows the probability density functions 
(PDFs) of the local shape factors H v (x,t) of the z- 
profiles of the instantaneous local velocity, obtained at 
x/D = 0.1, 0.2, 0.3, and 0.4 within region I. The dashed 
vertical line in the figure denotes the Prandtl-Blasius 
value for comparison. It is seen that the distributions 
are exactly peaked at the Prandtl-Blasius value, except 
at the position x/D = 0.4, where the peak is slightly 
off. This illustrates that most of the time the instan- 



taneous local velocity profiles in region I are indeed of 
Prandtl-Blasius type. Note that the position x/D = 0.4 
is close to the boundary x a between the regions I and II 
and hence the peak's slight deviation from the Prandtl- 
Blasius value at this position is likely to be caused by the 
competition between the corner roll and the LSC. 

Figure injb) shows the PDFs of H v (x,t) obtained at 
x/D = 0.5, 0.6, and 0.7 within region II. One sees that 
the peak positions of the distributions move closer to 
H^ B when proceeding from position x/D = 0.5, where 
the LSC is still slighly tilted downwards flowing, to the 
more plate parallel flow position of the LSC at x/D = 0.7. 
Specifically, at position x/D = 0.7, the PDF is exactly 
peaked at the Prandtl-Blasius value. This illustrates that 
the z-profilcs of the instantaneous local velocity are be- 
coming more Prandtl-Blasius along the evolution of the 
LSC. Referring to Fig. [T]this may be understood as fol- 
lows. The Prandtl-Blasius BL starts to develop in this 
region from the boundary between the corner roll and the 
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FIG. 11. (color online). The horizontal (independence of the 
deviations of the local velocity profile shape factors from the 
Prandtl-Blasius value for the laboratory frame, 8H v (x) (open 
triangles) and the time averaged dynamical frame, 5Hy(x) 
(solid circles). The vertical solid lines mark the boundaries 
between the three regions I, II, III. 



LSC near x/D = 0.4, and as one moves downstream the 
LSC becomes stronger and steadier so that it produces a 
more Prandtl-Blasius-like laminar layer. - The PDFs of 
H v (x, t) measured at x/D = 0.8 and 0.9 within region III 
are plotted in Fig. |^c) . Again, the two distributions are 
peaked close to H^ B , indicating that the instantaneous 
local velocity profiles obtained at these positions are of 
Prandtl-Blasius type for most of the time. 

We now consider the z-profiles of the instantaneous 
local horizontal velocity in direct comparison with the 
profiles of the times averaged fields seen from the lab- 
oratory system. Figure QJj] shows this direct comparison 
among the velocity profiles obtained at nine different hor- 
izontal positions xc the dynamical frame based local in- 
stantaneous horizontal velocity u*(x,z*), the laboratory 
frame based time-averaged velocity profile u(x, z), and 
the Prandtl-Blasius kinematic BL profile. Overall, ob- 
viously the u*(x,z*) profiles obtained in the dynamical 
frames match the Prandtl-Blasius profile well. This can 
be understood from our results in Fig. |H]that most of the 
time the instantaneous velocity profiles are of Prandtl- 
Blasius type and hence averaging all the rescaled profiles 
in the dynamical BL frames would naturally yield a pro- 
file of Prandtl-Blasius type. 

Figures [TOT a) to (d) show the profiles obtained in re- 
gion I. Both \u(x, z)\ and u*(x, z*) approximately match 
the Prandtl-Blasius profile for the range z* < 2, i.e., in 
the proper BL range. This suggests that the plume flow 
and temporal dynamics of the BLs do not play a key role 
in this region, which may be attributed to the strong 
mixing between the corner roll and the LSC. Note that 
after reaching the maximum values, both |u(x,z)| and 
u*(x, z*) decrease towards the bulk of the closed con- 
vection cell, while the Prandtl-Blasius profile keeps un- 



changed because it describes the situation of an asymp- 
totically constant, non-zero flow velocity. 

FiguresfTOTe) to (g) show the velocity profiles measured 
in region II. One observes that the time-averaged profiles, 
\u(x, z)\, obtained in the laboratory frame is much lower 
than the Prandtl-Blasius profile in the region around the 
kinematic BL thickness. As discussed in our previous 
papers^, a simple average of velocities at a fixed height 
z in the laboratory frame will sample a mixed dynamics, 
one pertaining to the BL range and the other one per- 
taining to the bulk, owing to the fluctuations of the BL 
thickness, and thus will distort the shapes of the profiles 
from that of Prandtl-Blasius. In contrast, the u*(x,z*) 
profiles measured in the instantaneous local dynamical 
frame agree well with the Prandtl-Blasius velocity pro- 
file, suggesting that the dynamical BL rescaling method 
can effectively disentangle the mixed dynamics of the BLs 
and the bulk, as all profiles arc expressed in the intrinsic 
BL-lcngth scale. 

Figures fTOTh) and (i) show the velocity profiles ob- 
tained in region III. Because of the small negative part of 
the profiles near the plate (see Fig. [He)), we better plot 
u*(x, z*) instead of \u(x, z)\ in the figures for comparison. 
One observes that u*[x, z*) averaged in the dynamical 
BL frames also here agree well with the Prandtl-Blasius 
z-profile, despite the existence of very complicated corner 
rolls in this region. 

The quantitative deviations of the velocity profiles 
shown in Fig. [TU] from the Prandtl-Blasius profile in 
terms of the shape factors are plotted in Fig. [TTJ 
When time-averaging in the laboratory frame, the veloc- 
ity shape- factor deviations 5H v (x) are closed to near 
the cell sidewall (x/D = 0.1, 0.2, and 0.3), but far away 
from at the other positions. In contrast, the shape- 
factor deviations 5H*(x) of the velocity profiles in the 
dynamical time-dependent local frames obviously are all 
much closer to zero, except that at x/D = 0.1, which is a 
bit of over-corrected. These quantitative results further 
indicate that our dynamical BL rescaling method works 
well for nearly all horizontal positions. And they con- 
firm the interesting fact that the Prandtl-Blasius laminar 
BL theory in RB flow is well justified despite the signifi- 
cant time dependence due to a vivid plume dynamics for 
Rayleigh numbers Ra below the ultimate state. 



IV. TEMPERATURE PROFILES NEAR THE PLATE 

We now turn to the corresponding analysis of the tem- 
perature field. In Fig. Q2]we compare the time-averaged 
thermal BL z-profiles obtained in the laboratory frame 
at nine different horizontal positions with the Prandtl- 
Blasius profile. The vertical distance z is normalized 
by the respective local thermal BL thickness Xth{x) and 
the temperature gradient of the Prandtl-Blasius profile 
is matched to those of the thermal BL profiles. Within 
region I (Fig. [T27a)). the profiles approach the bulk tem- 
perature faster, i.e., they run more quickly to the mean 
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FIG. 12. (color online), (a) The z-profiles of the time- 
averaged temperature field 0(x, z) as functions of the normal- 
ized distance z/Xth(x) using the lab- frame local thicknesses, 
obtained in (a) region I at x/D — 0.1, 0.2, 0.3, and 0.4; in 
(b) region II at x/D = 0.5, 0.6, and 0.7; in (c) region III at 
x/D = 0.8 and 0.9. The dashed lines indicate the Prandtl- 
Blasius thermal profile for comparison. 



bulk temperature = 1 when away from the sidewall. 
This is because most of the hot plumes rise upwards along 
the sidewall and hence the influences of thermal plumes 
on the thermal BL profiles become weaker when away 
from the sidewall. Figure IT27 b) shows the thermal BL 
profiles obtained at three positions within region II. One 
can notice that the agreements with the Prandtl-Blasius 
profile become worse as one moves from the still tilted 
down flow range of the LSC at the position x/D = 0.5 
to its more plate parallel flow at x/D = 0.7. This trend 




FIG. 13. (color online). PDFs of the shape factors of the 
rescaled instantaneous local temperature profiles (profiles sat- 
urate at first maximum, see text for explanation) obtained at 
x/D — 0.1, 0.2, 0.9. The dashed line marks the shape 
factor of the Prandtl-Blasius thermal profile for comparison. 



is different from what was observed for the kinematic BL 
profiles shown in Figs. [8jb) . - Finally, the temperature 
profiles obtained in region III are shown in Fig. [T27c). 
Both profiles are significantly lower than the Prandtl- 
Blasius profile for z/X t h(x) > 0.5, which we also attribute 
to the rising hot plumes in this region. 

To calculate the shape factors of the instantaneous lo- 
cal temperature profiles, Hth(x, t), we note that the emis- 
sions of plumes from the thermal BLs would lead to a 
much slower approach of the temperature profiles to the 
asymptotic value, and thus lead to a much lower value of 
the shape factor, because the temperature adjacent to the 
BL is not able to immediately relax back to the bulk value 
when a plume is detaching from the thermal BL. This 
holds for both the time-averaged as well as the instanta- 
neous profiles, especially for those obtained in the regions 
near the cell sidewall. For example, see the instantaneous 
temperature profile obtained at x/D = 0.2 (red circles) 
in Fig. |U here the temperature at the edge of the BL 
only reaches about 90% of the asymptotic value. There- 
fore, when calculating H t h{x,t) we take the approach of 
Zhou et at£: the first maximum temperature near the 
plate is defined as the asymptotic value of the z-profile 
and is used to normalize the profile. 

Figure [T3l shows PDFs of the instantaneous local ther- 
mal shape factors H t h(x,t). In the figure, Hf h B = 2.61 
is also plotted as the vertical dashed line for comparison. 
Unlike the case of H v (x,t) in Fig. HI the distributions of 
the Hth(x, t) are nearly independent of the horizontal po- 
sition, i.e., they all collapse on top of each other, except 
those obtained very near to the two sidewalk (x/D = 0.1 
and 0.9), which are shifted a bit to the left. In addition, 
the distributions are all peaked close to Hf h B , indicating 
that most of the time the instantaneous local tempera- 
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FIG. 14. (color online). Comparison between temperature ^-profiles obtained at x/Z) = 0.1, 0.2, 0.9 near the bottom plate: 
dynamical Q*{x,z^ h ) (red solid lines), laboratory Q(x,z) (blue solid lines), and the Prandtl-Blasius laminar thermal profile 
(black dashed lines), see also figure [151 



ture profiles over most part of the plate are of Prandtl- 
Blasius type. 

Direct comparisons of the z-profiles of the time- 
averaged local temperature Q(x,z) and of the dy- 
namically rescalcd field Q*(x, z^ h ) with the theoreti- 
cal Prandtl-Blasius temperature profile at nine differ- 
ent horizontal positions are plotted in Fig. Q31 Overall, 
the mean profiles obtained in the dynamical BL frame, 
0*(x, z^ h ), are much closer to the Prandtl-Blasius profile 
than the laboratory frame based profiles, Q(x,z), for all 
horizontal positions x, especially for those obtained near 
the plate's center (x/D = 0.4, 0.5, and 0.6), which match 
the Prandtl-Blasius profile exactly. 

One noticeable feature of Fig. [I4]is that the mean tem- 
peratures obtained near the cell's sidewall (x/D ^ 0.3 or 
x/D ^ 0.7) are much lower than the mean bulk tem- 
perature = 1 even at positions far away from the 
proper BL range, for both Q(x, z) and Q*(x, z^ h ). As 
discussed above, we attribute this to the emissions of 
thermal plumes, which would lead to a much lower value 
of the shape factor. Therefore comparison of the shape 



factors of such profiles with the Prandtl-Blasius value is 
somewhat meaningless. If, however, we choose the tem- 
perature at some position outside of the thermal BL, such 
as z/X t h(x) = 3 or z^ h = 3, as the asymptotic value for 
these positions (rather than the global bulk value) and 
use the asymptotic temperature to normalize the profiles, 
the re-scaled mean profiles look quite different. Indeed, 
following this procedure, we can eliminate the influences 
of plume emissions on the thermal BL profiles. The ob- 
tained mean temperature profiles in the next figure are 
much closer to the Prandtl-Blasius type. 

Figurc[T5]shows the direct comparison between the var- 
ious temperature profiles, rescaled in the described way: 
the dynamical local profiles Q*(x, Zf h )/Q*(x, z% h = 3), 
the laboratory profiles Q(x,z)/Q(x,z = 3\th), and the 
Prandtl-Blasius laminar thermal BL profile. Again we 
find a significant preference of the dynamical frame based 
profiles: Around the thermal BL thickness the laboratory 
frame based time averaged local profiles Q(x, z) /Q(x, z = 
3Ath) are all much lower than the Prandtl-Blasius profile, 
while the dynamically rescalcd instantaneous local pro- 
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FIG. 15. (color online). Comparison between specifically (see text) normalized vertical temperature profiles obtained at x/D = 
0.1, 0.2, 0.9 near the bottom plate: dynamical Q*(x, z* h )/Q*(x, z% h = 3) (red solid lines), laboratory Q(x,z)/Q(x,z = 3Xth) 
(blue solid lines), and the Prandtl-Blasius laminar thermal profile (black dashed lines), see also figure [TH 



files Q*(x, Zt h )/@*(x, z^ h = 3) match the Prandtl-Blasius 
profile much better. The shape-factor deviations of these 
rescaled profiles are shown in Fig. [16] The shape-factor 
deviations SH t h{x) for the laboratory frame profiles are 
definitely smaller than zero. In contrast, the dynamical 
frame based deviations 5H^ h (x) are much closer to zero, 
suggesting that our dynamical BL rcscaling method can 
indeed capture the BL properties efficiently 

V. CONCLUSIONS 

In conclusion, we have made a systematic study of the 
horizontal (a;-) dependence of the shapes of the z-profiles 
of the kinematic and thermal BLs in turbulent RB con- 
vection using 2D numerical data. We have extended our 
previous studies, which were restricted to the plate's cen- 
ter, to all horizontal positions along the bottom (or top) 
plate. The major findings can be summarized as follows: 

1. In situations where the traditional methods based 
on the time-averaged horizontal velocity profiles are 



no longer capable of producing physically mean- 
ingful BL thicknesses, the time-averaged instanta- 
neous BL thicknesses provide well-defined length 
scales for both the kinematic and thermal BLs. 
Such situations can arise, for instance, from the 
competition between the LSC in the center region 
and the secondary rolls near the corners. 

2. When the instantaneous local velocity and temper- 
ature values are rescaled by their respective instan- 
taneous local BL thicknesses, it is found that the 
Prandtl-Blasius profiles hold in an instantaneous 
sense most of the time. 

3. For most parts of the horizontal bottom (or top) 
plate both the local velocity and temperature pro- 
files match the classical laminar Prandtl-Blasius BL 
profiles well, if they are re-sampled in the respective 
dynamically rescaled frames, which fluctuate with 
the instantaneous local kinematic and thermal BL 
thicknesses. 
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FIG. 16. (color online). The horizontal (x-)dependence of the 
deviations of the local temperature profile shape factors from 
the Prandtl-Blasius value for the laboratory frame, 8H t h(x) 
(open triangles) and the time averaged dynamical frame, 
8Ht h (x) (solid circles). Here, both SH t h{x) and 8H^ h (x) are 
calculated based on the profiles in figure[T5] The vertical solid 
lines mark the boundaries at x a and Xb between the three re- 
gions I, II, III. 
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